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Formal solutions to the mode equations for both spherically symmetric black holes and 
Bose-Einstein condensate acoustic black holes are obtained by writing the spatial part of the 
mode equation as a linear Volterra integral equation of the second kind. The solutions work 
for a massless minimally coupled scalar field in the s-wave or zero angular momentum sector 
for a spherically symmetric black hole and in the longitudinal sector of a ID Bose-Einstein 
condensate acoustic black hole. These solutions are used to obtain in a rigorous way analytic 
expressions for the scattering coefficients and gray-body factors in the zero frequency limit. 
They are also used to study the infrared behaviors of the symmetric two-point function and 
two functions derived from it: the point-split stress-energy tensor for the massless minimally 
coupled scalar field in Schwarzschild-de Sitter spacetime and the density-density correlation 
function for a Bose-Einstein condensate acoustic black hole. 
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I. INTRODUCTION 

A major prediction of quantum field theory in curved space is that when a black hole, 
BH, forms from the gravitational collapse of a star, particle production occurs and the 
BH evaporates [lj. If backreaction effects are neglected then, after a transient which 
depends on the details of the collapse, the emission rapidly becomes stationary, is nearly 
thermal, and depends only on the geometry exterior to the star and to the event horizon 
of the BH. The temperature of the radiation is proportional to the surface gravity k of 
the event horizon. The fact that this radiation depends only on the parameters (conserved 
charges) characterizing the spacetime exterior to the BH is a manifestation of the BH 
no-hair theorem (see for example 0). 

The spectrum of the radiated particles is not exactly thermal because the modes of the 
quantum fields propagating from the horizon region to infinity experience the presence of 
an effective potential that causes partial reflection back to the horizon. So, the number of 
particles which propagate to infinit}0for each mode (in units where ks = fr = c = 1) is 

jvC) = r( ‘P> 


27 TUJ 

e — 1 


( 1 . 1 ) 


Here k is the surface gravity of the event horizon, oj is the frequency of the mode, and 
i denotes the other quantum numbers needed to specify the mode. The factor T^i(o;) is 
the so called ‘gray-body factor’ which accounts for the backscattering mentioned above. 
In the literature rW(u;) is usually computed analytically by the method of asymptotic 
matching [3H3- Since the effective potential vanishes at the horizon and at infinity (or 
the cosmological horizon) it is possible to write the solutions to the mode equation at 
these locations in terms of plane waves. These are then matched in the intermediate 
region. Particular interest lies in the low frequency behavior of r^(o;) because the infrared 
divergence coming from the Planckian factor in (1.1) can be avoided if r®(cu) vanishes in 


1 In the case of a black hole immersed in de Sitter space this is the number of particles that propagate to 


the cosmological horizon. See Sec. VB for more details. 
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the oj —> 0 limit. If rW(cu) approaches a nonzero constant in this limit then an infinite 
number of ‘soft’ particles will be emitted. 

For asymptotically flat spherically symmetric BHs one has for the i = 0 mode, which is 
the dominant one for BH radiation, a universal behavior at low frequency r^ = °) (w) ~ Ahu 2 
where Ah is the area of the event horizon Pi- However for BHs in de Sitter space one 
has r^ =0 )(u;) ~ const, in the low frequency limit for the massless minimally coupled scalar 
field |10l TT1 j^] A nonvanishing rh)(cu) in the cu —> 0 limit was also found recently [I3J for 
a class of Bose-Einstein condensate, BEC, acoustic BHs which have effectively have one 
spatial dimension. It is important to stress that the infrared divergences in the number 
of created particles found for the Schwarzschild de Sitter and BEC cases are weak enough 
that there is no infrared divergence in the combined energies of these particles. 

A possible problem with the asymptotic matching method described above and used 
in pain is that the near-horizon and zero frequency limits do not commute for solutions 
to the mode equation. The same is true for the zero frequency and infinite distance limits 
(or in the Schwarzschild-de Sitter case, the zero frequency limit and the limit in which the 
cosmological horizon is approached). 

In this paper, we provide rigorous derivations of the zero frequency limits for the scat¬ 
tering coefficients and gray-body factors for the mode functions of the massless minimally 
coupled scalar field in spacetimes containing spherically symmetric black holes and for BEC 
acoustic black holes. For spherically symmetric black holes, we restrict our attention to 
modes with zero angular momentum, i.e. the s-wave sector. For BEC acoustic black holes 
we assume that the system can be approximated as having just one space dimension and 
thus work in the longitudinal sector. We also only consider cases in which the velocity 
v of the condensate is constant and the sound speed c varies as a function of position. 
Qualitatively similar results can be obtained for more general cases in which both c and v 
vary with position. 

2 For nonminimal coupling to the scalar curvature or i ^ 0 one recovers even for these BHs a vanishing 
gray-body factor (see for example [T2J and US). 
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To compute the scattering coefficients and gray-body factors we have developed a gen¬ 
eral analytic method, based on a linear Volterra integral equation of the second kind [T5], 
for solving the mode equation. In all cases, we compute the exact (to leading order in 
u>) scattering coefficients in the exterior of the black hole and relate them to the horizon 
boundary values of the appropriate exact solutions (and their first derivatives) of the zero 
frequency mode equation. 

We use the results to investigate infrared divergences in the symmetric two-point func¬ 
tion for the massless minimally coupled scalar field in the Unruh state [T6j in both the 
BH and BEC cases. To our knowledge the only prior work in this area is in (14] where 
it was shown that an infrared divergence occurs for the two-point function in the BEC 
casej^] We find that the two-point function is infrared finite if the gray body factor van¬ 
ishes in the zero frequency limit. However, it is infrared divergent if the gray body factor 
does not vanish in the zero frequency limit. From the relationship between the two-point 
function and the point-split stress-energy tensor for a spherically symmetric black hole or 
the density-density correlation function for a BEC acoustic black hole, one would guess 
that infrared divergences in the two-point function would lead to infrared divergences in 
these quantities as well. We show that this is not the case and that there are no infrared 
divergences in either the point-split stress-energy tensor for the massless minimally coupled 
scalar field in Schwarzschild-de Sitter spacetime or the density-density correlation function 
for the class of BEC acoustic black holes that we consider. 

In Sec. [TT] we find two formal linearly independent solutions to the mode equation in 
terms of Volterra integral equations of the second kind. We then find a second set of 
formal solutions to the spatial part of the mode equation which are useful near the event 


horizon. In Sec. Ill we use these solutions to derive formal expressions for the scattering 


coefficients and gray-body factors at all frequencies. In Sec. IV it is shown that in the low 
frequency limit the scattering coefficients and gray-body factor can be written in terms of 


3 There were also results announced in [14] for the two-point functions for spherically symmetric black 
holes, for the point-split stress-energy tensor, and the density-density correlation function, but the proofs 
of these results are given here. 
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exact solutions to the zero frequency mode equation which are evaluated in the limit that 
the event horizon is approached. The condition which determines whether the gray-body 
factor vanishes or approaches a nonzero constant in the zero frequency limit is given. In 
Sec. [V] specific expressions for the low frequency limits of the scattering coefficients and 
gray-body factor are derived for Schwarzschild, Reissner-Nordstrom, and Schwarzschild- 
de Sitter spacetimes as well as for the class of BEC acoustic BHs we are considering. In 


Sec. VI the question of infrared divergences in the two-point function, the point-split stress- 
energy tensor, and the density-density correlation function are investigated. A summary 


and discussion of our results is given in Sec. VII In Appendix A some of the details of the 


derivations in Sec. IV are given. In Appendix B exact (in w) near-horizon solutions to the 
mode equation are constructed. 


II. VOLTERRA INTEGRAL EQUATION 

The Klein-Gordon equation for a massless minimally coupled scalar field 

-J=d ll (V-gd' , ^)=o (2.1) 

arises both for a quantum field in curved space and as a description of the quantized phase 
fluctuations (phonons) of a BEC in the hydrodynamic limit (details are given, e.g., in jl7j). 
For both an analog black hole and a gravitational black hole, the class of metrics we shall 
be interested in can be written in the form 

dr 2 

ds 2 = —f(x)dt 2 + + C\(x)dx]_ , (2.2) 

I \ x ) 

where / = 0 (and x = 0) at the horizon, and we have explicitly separated the (1 + 1) and 
‘transverse’ parts. For a spherically symmetric black hole 

Cl(x)dx 2 ± = (x + r H ) 2 dn 2 , (2.3) 

with rii the radius of the event horizon. The usual radial coordinate is r = x + ru and 
we use units such that h = c = G = kg = 1. For a BEC it can be shown that the original 







7 


acoustic metric can be mapped, via a coordinate transformation, to the form (2.2) with 

(2.4) 


n 


Cl(x)dx 2 ± = —( dy 2 + dz 2 ) , 
me 

where n is the number density of particles in the condensate and m is the mass of an 
individual particle. Here we consider for simplicity only models in which the BEC moves 
with constant velocity v < 0 in the x direction, n is a constant, and the sound speed c varies 
with x (see e. g. but our analysis can be easily extended to models in which c, v and n 
vary with x. In what follows we ignore the irrelevant constants n and m which is equivalent 
to choosing units such that n/m = 1. Further in the gravitational case we consider only 
the spherically symmetric l = 0 modes and in the BEC case only the longitudinal ones. 
Thus in both cases cj) = <j>(t,x). For stationary solutions 

4>(x,t) = (j) u (x)e~ lut , 


and Eq. (2.1) becomes 


1 _d_ 

C 2 , (x) dx 


Ci(x)/(x)4 


OJ 


4>u(x) + = 0 . 

/(*) 


(2.5) 


( 2 . 6 ) 


A useful alternative form for the mode equation can be derived by making the changes 
of variable 

X(t, x ) 


4)(t, x) = 


X = 


C±(x) ’ 

r dx 


f(x) 


(2.7a) 

(2.7b) 


where x* is the tortoise coordinate. To determine the behavior of x* near the horizon note 


that for the cases we consider it is possible to expand the metric coefficient f(x) in (2.2) 
in powers of x such that 


f ( x ) = & 

The definition of the surface gravity is 


i— 1 


k = 1 df(x) 


2 dx 


( 2 . 8 ) 


(2.9) 


x=0 
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thus near the horizon 


/ = 2kx + 0(x 2 ) . 


( 2 . 10 ) 


Substituting into (2.7b) it is easy to see that near the horizon x* —> — oo and more precisely 
that 


x = xoe 


2 KX* 


+ 0(e 


4kx* ) 


( 2 . 11 ) 


with xq an arbitrary positive constant. Far from the horizon, if / —> const, as x* —> +oo 
then x* ~ x. If instead there is a cosmological horizon then x* —> +oo at the cosmological 
horizon. 


Substituting (2.7) into (2.1) gives 


{-dt + d 2 x , - H e ff] x(t, x) = 0 . 


( 2 . 12 ) 


Here we have introduced the effective potential 


„ 1 <?C ± 1 , d (,dC L \ 


C± d(x*) 2 C ± J dx 

For stationary solutions x(x,f) = Eq. (2.12) takes the Schrodinger-like form 

d 2 


(2.13) 


+ W — V e g 


Xu = o 


(2.14) 


d(x*) 2 

For the cases we consider V e ^ vanishes at x* = ±oo. Thus the general solution to 


Eq. (2.14) will be a linear combination of two solutions Xu an d Xu with the asymptotic 


behaviors 


Xa 


COSUJX 


Xu 


smwi 


as 


oo . 


(2.15) 


Formal solutions for Xu an d Xu can obtained in terms of the following linear Volterra 
integral equations of the second kind: 

^ /*°° 


Xu( x *) = COSU1X* 


S / * \ • * 

Xu\x ) = sin wx 


u J x . 

i r°° 
u .L* 


sin[w(x* - y*)}V eff (y*)x C u(y*) , 
sin[w(x* ~ y*)}V eff (y*)x S u(y*) ■ 


(2.16a) 

(2.16b) 
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For all x* these solutions can also be written in the general form 


Xu( x *) = A{oj, x *) cosujx* + B(u, x*) sincux* , 

XLi. x *) = C{u,x*) cosux* + D(u, x*) sin lux* , (2-17) 


with 


A[u,x*) = 1 + - [ dy* sin uy* V eS (y*)xl(y*) , 
a; J x * 

i r°° 

B(oj,x*) = -~ dy* cos uy* V e s(y*)xZ{y*) , 

W Jx* 

C(u,x*) = ~ [ dy* sinuy* V e s(y*)xtj{y*) , 

1 Z* 00 

D(oj,x*) = 1-/ dy* coswy* F eff (y )x® (y ) • 

^ JX* 

We note that 


dx* 

dxS, 

dx* 


—^4(o;, x*)w sinwx* + -B(w,x*)u;coswx* , 
—C(a;,x*)wsinu;x* + .D(w, x*)w coswx* , 


and that the Wronskian 


W = 


* dxuj 
dx* 


dxt, 

dx* 


Xu , 


(2.18) 


(2.19) 


( 2 . 20 ) 


is a constant. For Xu = Xu + it can easily be seen by evaluating (2.20) in the large x* 


limit, that W = 2 ioj and thus that 


1 

ui 


Xu 


dxi 

dx* 


X S u 


dxi 

dx* 


= 1 


Using Eqs. (2.17) and (2.19) one finds 


( 2 . 21 ) 


A(ui,x*)D(ui,x*) — B(oj,x*)C(u:,x*) = 1 . 


( 2 . 22 ) 
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III. SCATTERING COEFFICIENTS AND GREY-BODY FACTORS 


In the horizon limit, x* -A —oo, the solutions (2.17) have the asymptotic behaviors 


Xu 


Xlx 


Aj cos lux* + Bu sin lux* 
C u cos lux* + Du sin lux* 


(3.1) 


with 


1 r°° 

Au = A(lu, -oo) = 1 + - / dy* sin buy* V eS (y )xZ(y ) , 

Bu = B(u, -oo) = - / dy* cosc uy* V e g(y*)Xu(y*) ■ 

l r°° 

Cu = C(u, -oo) = - dy* sin toy* V eS (y )xL{y ) , 

^ J—oo 

1 z -00 

Aj = -oo) = 1-/ dy* cos wy* V e g(y )%* (y ) . 

^ J — OO 


(3.2) 


We shall now relate the coefficients (3.2) to the scattering coefficients in the exterior of 


the black hole (the / > 0 region in (2.2)) [17j . Because of the independence of Hawking 
radiation from the details of the collapse producing the BH, when discussing this emission 
one can consider just the spacetime exterior to the star, mimicking the effects of the collapse 
by giving appropriate boundary conditions on the past horizon of the analytically extended 


exterior manifold. A complete basis for the solutions of the mode equation (2.12) is formed 
by two sets of modes, xi an d xh, shown in Figs. 1, 2 |18j . The modes xi originate at past 


null infinity I (x* = — oo) and because of the potential term in Eq. (2.14) are partially 


transmitted towards the future horizon H + (x* = +oo) and partially reflected to future 
null infinity I + ( x* = +oo). The modes xh originate on the past horizon H~ ( x* = — oo). 
They are partially transmitted to future null infinity I + , and partially reflected back to 
the future horizon H + . More specifically 


Xi(t,x ) 

= A e -<w x / (l) , 

\JLU 

(3.3a) 

XH(t,X ) 

= A=e-‘“‘ x »(x) . 

V uj 

(3.3b) 
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FIG. 1. Modes Xh starting at past null infinity (/ ), transmitted to the future horizon ( H + ) and 
reflected to future null infinity ( I + ). 



FIG. 2. Modes xh, starting at the past horizon (H ), transmitted to future null infinity ( I + ) and 
reflected to the future horizon (H + ). 


with xL and Xu solutions to Eq. (2.14). For a ID BEC acoustic black hole N = and 


for a spherically symmetric black hole N = 


It is useful to write \L an d Xu dr terms of solutions to the mode equation (2.14) which 
correspond to outgoing (e lux ) and incoming (e~ luJX *) waves at x* = oo. These solutions 
are 


X? = Xu + ixl 


xT = x c u 


ixt 


(3.4) 
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They have both right and left moving parts near the horizon (x* —>■ — oo) so that 

ITT iux* , 7-1 —iux* 

hi r e + , 




Xz 


Eie %ux * + F ie ~ 


(3.5) 


where, using (3.1) 


E r = — [via; + Du — i(Bu — C(j)] , 
Er = ^ [j4cj — Du + i(-Ba; + C^)] , 
Ei = — [Au — Du — i(Bu + Co;)] , 
Ei — — [Au + Du + i{B u — Cu)} ■ 


(3.6) 


To construct xL we consider a mode corresponding to an initial unit amplitude left- 
moving wave at x* —> +oo which gets partially reflected back to x* —> +oo and partially 
transmitted to the horizon (x* —> —oo). Then 


xL = xT + R ix 


oo ^ ^—iux 


+ Rie 1 ' 


as x 


+oo , 


Near the horizon this solution has the behavior 


Xu 


Tje: 


(3.7) 


(3.8) 


We find, from (3.5) and (3.6) that 


u E[ (Au — Du — i(Bu T Cu)) 

RI = -~FT = —- 


E r 


Tj = F X - 


EiF r 


E r 


1 

E r 


(Au T Du — i(Bu — Cu)) 

2 

(Au + Du- i(Bu - Cu)) ’ 


(3.9) 


where in the last equation we have used the Wronskian condition (2.22). 

To construct Xu we s t ar t with a mode corresponding to a unit amplitude right-moving 
wave which originates at the horizon (x* —> — oo) and is partially reflected back to x* —> —oo 
and partially transmitted to x* —> +oo. Thus 


Xu = T HXr 


ThC 


as x 


Too . 


(3.10) 
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Near the horizon this solution has the behavior 




+ -R/f e 


(3.11) 


Thus we find, from (3.5) and (3.6) that 

p _ Fr _ (Aj — D u + + C^)) 

H ~K~ (A, + Rla, - Z(S U - a)) ’ 
1 9 
Th = -e = 


{Auj + .Du — i(Buj — Cuj)) 


Finally, the gray-body factor is defined as 


(3.12) 


r=| t h \‘ 


(3.13) 


and represents the probability that a mode which originates on the past horizon H and 


has unit norm reaches I + ( x* = +oo). It appears in Eq. (1.1) where it modifies the exact 
black-body spectrum for the number of particles which reach I + at a given frequency uj. 


IV. LOW-FREQUENCY LIMIT 


In this section we consider the low-frequency limits of the expressions derived in the 


previous subsections. Therefore we need to analyze, the behaviors of the integrals (3.2) for 
small values of u;. The main difficulty with evaluating these integrals in the limit uj —> 0 is 
that this limit does not commute with either of the limits x* —>• Too and it is not possible 
to compute the integrals analytically for arbitrary values of u. However, one can break the 
integrals into three parts so that 


f 


dy* f A dy* + ^ dy* + 
J— oo J— A 


dy* 


(4.1) 


and take the limits uj —> 0 and A —> oo in such a way that uA <C 1. Then expansions 
can be used near the horizon and at large y* for the potential V e s in the first and third 
integrals which make it possible to compute these integrals analytically. It is shown in 
Appendix A that when this is done these integrals vanish for all values of uj including 
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d = 0 provided V e g approaches zero fast enough at the horizon and at infinity. It is also 


shown in Appendix A that for the middle integral we can use Eq. (2.14) to obtain the 
relation V e g{y*)Xu{y*) = [dy* + d 2 ]xw(y*), integrate twice by parts, and obtain 


A u = 1 + lim 

A—>-oo 


B u = — lim 

A—>-oo 


Coj = lim 
A—>-oo 


sin(t<;y*) d 
d dy* 

cos (dy*) d 


-cos (dy*) ) xliv*) 
i A 


A 

J -A 


d dy 
sin(wy*) d 


+ s in(wy*) xZ(y*) 


J -A 


D u = 1 — lim 

A—)-oo 


d dy 

cos(c dy*) d 
d dy* 


- cos(dy*) xl(y*) 


A 


J -A 


+ sin (dy*) xZ(y*) 


-i A 


(4.2) 


J -A 


Using the asymptotic behaviors (2.15) it is easy to show that in the limit A —> oo 
these coefficients only depend on the parts of the above expressions which are evaluated at 
y* = —A. Thus 

sin(u;y*) d 


A u = lim 

A —>oc 


= lim 

A —>OC 


Cu = lim 
A—>-oo 


D u = lim 

A—>-00 


cos {uy*) — 

cos (uiy*) d 
id dy* 

cos (uiy*) — 


d dy* 


xKv* 


+ S\n(idy*) xZ(y*) 


sin(u;y*) d 


d dy* 




cos(dy*) d 


+ sin (uy*) xt(v*) 


y *=-A 

y *=-A 

y *=-A 

y *=-A 


(4.3) 


d dy* 

Recalling that we take the id —> 0 and A —> oo limits in such a way that id A < 1, we define 
the following constants: 


= Aq = lim 

u ->0 

A —>oc 


sin(cm/*) d . 
cos (dy )-— Xu\y 


AS = lim (dBj) = lim 

OJ — oj —^0 

A—>00 L 

‘S? = lim(u; _1 C w ) = lim 

OJ — oj —>-0 

A—»oo L 


id dy 
cos (dy*) d 


d dy* 


y *=-A 
+ sin (dy*) ) dxZ(y* 


J y *=-A 


cos(u;y*) — 


sin (uy*) d \ xl(y*) 


d dy* 


d 


y *=-A 


= Dq = lim 

oj— 

A—>-oo L 


cos(o;y*) d 
d dy* 


+ sin (dy*) xt(y*) 


- y *=-A 


(4.4) 
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Using the Wronskian condition (2.21) it is easy to show that 


- SBV = 1 . 

We fined that for small cu 

= stf + 0(cu 2 ) , 
C3) 

Bu = ~ + 0(u) , 

CJ 

Cu: = u'io + O(lu^) , 
D u = ® + 0{ cu 2 ) . 


(4.5) 


(4.6) 


To evaluate the expressions in (4.4) we first note that since cuA -C 1 


SB 

9 


lim 

u^O 

A—>-oo 

lim 

A —>oo 


i-v'A)x c Jv') 


y *=-A 


dy* 


xZiv*^ 


y *=-A 


lim 

w->-0 

A—^oo L 


1 _ y*—\ AUy*) 

dy*) cu 


Jy*=-A 


lim 

u->0 

A—>-oo 


d xZ(y* 

dy* cu 


Jy*=-A 


(4.7) 


Next we show that both xZ( x *) an d *“u ^ and their first derivatives can be evaluated in 
the limit cu —> 0 with uj\x*\ <C 1 by identifying them with solutions to the mode equation 
when cu = 0. To do so we consider small values of cu and large values of x* with two 
constraints. One is that cu 2 » V e s(x*) which means that to leading order xZ = cos(c dx*) 
and xto = sin(cu.x*). The second is that 0 < uix* <C 1 which means that to leading order 
Xu = 1 and Xu = ux *. Examination of the large x* behaviors given in Appendix A for 
V eS in the cases we are interested in shows that it is not hard to find values of cu and 
x* that satisfy these two conditions. We want to match xZ an d Xu to solutions to the 


4 The relation for A u is proved in Appendix A. The proofs of the other relations are straight-forward 
generalizations of that proof. 
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mode equation (2.14) with uj = 0. We will denote the solution that xZ matches to as 

The situation for xt is more subtle because xZ ~ <xx* and thus vanishes in the limit 

u —> 0. However, ^ ~ x* does not vanish in this limit. Further this quantity obeys the 

same mode equation as x'Z and in the small a; limit this is just the zero frequency mode 

equation. Therefore we can match the large x* behavior of with the corresponding 

( 2 ) 

behavior of a solution to the zero frequency mode equation which we will denote by Xo ■ 
Thus, in the limit x* —> oo these solutions have the behaviors 


Xo\x*) -> 1 • 
Xo\ x *) -> x* 


(4.8a) 

(4.8b) 


To finish the matching we need a condition on the derivatives of these solutions. This can 


be obtained from the Wronskian condition (2.21) satisfied by xZ and xZ- The corresponding 
condition for Xq 1 ’ 2 ^ is 


, ( 2 ) j (!) 

v (l)^V _ (2) ^Xo_ = 

X ° dx* Xo dx* 


(4.9) 


This fixes the asymptotic behaviors of their derivatives. Then we can write for u —> 0 with 
a;lx*I <C 1 


V c (1) 

lim n Xw = Xo 

u —>-0 

i* Xuj (2) 

Imr — = Xo 

0 UJ 


For future reference we note that therefore 


,. oo (1) , . (2) 

hm x° =X 0 + lU} Xo 

'jJ —^0 

lim i xT = X (1) - ivxb 2) ■ 


(4.10) 


( 4 . 11 ) 
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f 1 2) _ 

Replacing xZ an d xt]/ u with the corresponding solutions Xq ’ ' in Eq. (4.7) gives 






hor 


hor 






- hor 


(4.12) 


- hor 


where the subscript “hor” means that the quantities inside the square brackets are to be 


evaluated on the horizon. Note that the Wronskian condition (4.9) implies that 


- W = 1 . 


(4.13) 


Substituting (4.3) into (3.9) and (3.12) and using (4.6) one can see that at low frequency 


the coefficient SB plays a crucial role. If SB ^ 0 then 


2iu 


Ri, Rh = —1 + 0(uj), T i = T h = — + 0(uj z ), 


(4.14) 


while if SB = 0 we have the very different behaviors 


Ri = 


+ 


+ 0(u>), Rh — 


+ 


+ 0(u>), T i = T h = 


sf + 3> 


+ 0(u) . (4.15) 


Similarly the gray-body factor (3.13) vanishes at low-frequency as 

r 


4w 2 

w 


(4.16) 


when SB ^ 0, while it approaches the constant value 

r = 4 


(4.17) 


{st + 

when SB = 0. 

From Eq. (4.12) it is clear that SB = 0 if /dx* vanishes at the horizon. The mode 


equation (2.14) can easily be solved exactly when uj = 0 by going back to Eq. (2.6), setting 
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uj = 0 and using the relation 4>q = xo /Cj_. The general solution is 

dx 


Xo{x) = clC_l(x) + b 
with a and b arbitrary constants. Then 


C±(x) J 3 


C\ (x)f{x) ’ 


dx o = J_ dC± b 

dx* ~ C± dx 1X0 + C± 


(4.18) 


(4.19) 


Since for the models we consider, / vanishes at the horizon and C± and its first derivative 
are nonzero constants there, it is clear that dx^ /dx* can only vanish at the horizon if 
6 = 0. However, Xq 1 ^ is specified by the condition ( |4.8a ) which can only be satisfied with 
6 = 0 if C_ i approaches a constant in the limit x* —> oo. As shown below this occurs for 
Schwarzschild-de Sitter spacetime and for a ID BEC acoustic black hole. It does not occur 
for any spherically symmetric black hole spacetime that is asymptotically flat. 

In cases where the solution Xo has 6 = 0, it is necessary that V e s changes sign at some 
finite value of x*. To see this note that for the cases we consider near the horizon 


C±(x) = C ± ( 0) + 


dC±(x) 


dx 


x + 0(x 2 


(4.20) 


:r=0 


Using (2.11) it is easy to see that near the horizon dC±/dx* and d 2 C±/d(x*) 2 have the 
same sign. Thus for C±_ to approach a constant value as x* —> oo it is necessary that 


d 2 C±/d(x*) 2 changes sign. Then Eq. (2.13) shows that V e Q must also change sign. 


V. SPECIFIC CASES 


In the following subsections we provide specific examples of the two different behaviors 
which can occur for the scattering coefficients and the gray-body factor in the low frequency 
limit. 


A. Schwarzschild and Reissner-Nordstrom spacetimes 

The first, obvious, application is Schwarzschild spacetime, which will allow us to com¬ 
pare the results of our general procedure with those in the literature |5] . Referring to the 
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FIG. 3. Effective potential in Schwarzschild spacetime (M = 1) plotted as a function of x = v — Th- 


generic metric (2.2), we have with r = x + th, f = 1 —F" = x +r H an< ^ ^-L = f = x + r#. 
Here rjj = 2 M is the horizon radius and M is the mass of the black hole. The tortoise 
coordinate (2.7b) is x* = f y- = x + j-r + fn In (x/rn)- 

For this case, Eq. ( |2. 14 ) is the mode equation for modes in the s-wave sector of a 
massless scalar field J^] As illustrated in Fig. 3, the effective potential (2.13) 


V e ff = (1 — —)-w = 


x th 


r ' r J (x + vh) a ’ 


(5.1) 


has the shape of a positive barrier, i.e. it starts from 0 at x = 0 (x* = — oo), reaches a 
maximum for x = \rn and then goes to 0 asymptotically as x,x* —> oo. The fact that 
it does not change sign is related to the fact that C± is unbounded in the limit x* —> oo. 


From the discussion near the end of Sec. IV this means that SB ^ 0 and thus the gray-body 
factor vanishes in the limit uj —> 0. 

To calculate the scattering coefficients and gray-body factor we need to find the solutions 


( 12 ) • ... - 

Xq to the zero frequency mode equation. The general solution to this equation (4.18) in 


5 For both Schwarzschild and Reissner-Nordstrom spacetimes the scalar curvature R is zero so the mode 
equation is independent of the coupling to the scalar curvature. 
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this case is 


( a \ a S x + r ^) i 
Xo = a(x + ru) + o -m 


rn 


x + r H 


The boundary condition (4.8a) for Xn 1 is satisfied if a = 0 and b = — 1 so that 


Xo 

(l) 0 + r H ) 

Xo = - ln 

r H 


x 


x + r H 


(5.2) 


(5.3) 


Thus Xo diverges as — at the horizon. 


rH 


( 2 ) . . 

For Xo the boundary condition (4.8b) fixes a = 1 but leaves b undetermined so that 

x 


( 2 ) , A h( x + rH h 

Xo = x + r H + b -ln 

r H 


x + rn 


(5.4) 

( 2 ) 


Note that the Wronskian condition (4.9) is satisfied for arbitrary values of b. Thus Xo 
rn + near the horizon. We show below that the ambiguity in the value of b does not 
affect the low-frequency behaviors of the scattering coefficients. 


(4.12). The result is 


easy to compute the coefficients &/, 8$, 

, <jf, 9 in 

srf = 1 , 

(5.5a) 

88=-— , 

(5.5b) 

r H 


& = r H -b , 

(5.5c) 

9= — ■ 

(5.5d) 


ro 


Since 88 ^ 0, one finds from (4.14) that 


Rj, Rh ~ —1 + 0{lo ) , Tj = Tfj ~ —2 irjjio + 0(w 2 ) . 


(5.6) 


To leading order the gray-body factor is 

T ~ 4r^u; 2 , (5.7) 

in agreement with the results in [5]. We see immediately the important role played by the 
gray-body factor, which prevents waves with frequency cj < from reaching infinity. 
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This analysis can be easily extended to Reissner-Nordstrom spacetimes, where the re¬ 
sults are qualitatively similar to those for Schwarzschild spacetime. In this case 


/= 1 - 


r+ 


1 - 


(5.8) 


x + / \ x + r +i 

with with Q the electric charge of the black hole, r# = r+ = M + \JM 2 — Q 2 the outer 
(event) horizon, and r_ = M — \JM 2 — Q 2 the inner (Cauchy) horizon. The tortoise 


coordinate (2.7b) is 


x* = x + r+ — (r+ + r_) ln(r+ — r_) + 


(r+ ~ r-) 


ln(x) — 


r+ — r_ 


ln(x + r+ — r_) . (5.9) 


The general solution (4.18) of the zero frequency mode equation is 

x T r + 


Xo = a(x + r+) + b 


r_|_ — r_ 


In 


x + r + — r_ 


(5.10) 


Imposing the conditions (4.8) one finds that 

x + r + 


Xo 1} = 


In 

r + — r_ \ x + r_|_ — r_ 


( 2 ) . x + r+ 

Xn = x + r + + o-In 

— r_ 




r — r_ 


(5.11) 


Substituting into (4.12) gives the same expressions for iz/, and ^ as in the 


Schwarzschild case except that the event horizon is at ru = r+ and r + ^ 2 M unless 
Q = 0. Similarly the expressions for the scattering coefficients and the gray-body factor 
are the same as for Schwarzschild to leading order. 

B. Schwarzschild-de Sitter spacetime 

We next consider a case which differs qualitatively from the previous ones, namely a 
black hole immersed in an expanding de Sitter universe and described by the Schwarzschild- 


de Sitter solution, with the metric coefficients in (2.2) given by / = 1---and 


C± = r. Here A is the cosmological constant. For 0 < A < 9A I 2 the space-time has 



















22 



FIG. 4. Effective potential in Schwarzschild-de Sitter spacetime (M = 1 and A = 3 x 10 3 ) plotted 
as a function of x = r — rjj- 


two horizons at r = ru and r = rc (with rn < rc) which are, respectively, the black 
hole and cosmological horizons. The equation / = 0 has an additional negative solution 
(r = —ro < 0) which is unphysical. We define the coordinate x as before so that x = r — rjj. 


We shall consider the mode equation (2.14) in the region between the two horizons, 
0 < x < rc — rjj, which is the region where t is timelike and x is spacelike. From 


Eq. (2.7b) it can be seen that the tortoise coordinate has the values x* = —oo at the event 
horizon ( x = 0) and x* = +oo at the cosmological horizon (x = rc — rjj). V e Q is plotted 
in Fig. [4} In addition to the barrier which was also present in the Schwarzschild case, there 
is a well outside the barrier in which V e s < 0. Since C± = rc at x* = oo, the gray-body 
factor does not vanish in the zero frequency limit. 


To compute the low-frequency scattering coefficients and gray-body factor in this space- 

( 1 , 2 ) 


time we use (4.18) and (4.8) to find the mode functions Xq . Since the term proportional 


to b in (4.18) diverges as In |(x + rjj — rc)/rc\, the boundary condition (4.8a) is satisfied 
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if a = — and b = 0, so that 

rc > 


(l) x + r H 

Xo =- • 

r c 

At the cosmological horizon the surface gravity is kc = h\f (x)\ x=rc - rH and 


(5.12) 


2 kc 


In 


I® - ( r c - r H ) | 


rc 


(5.13) 


Near the cosmological horizon, the term proportional to b in (4.18) has the leading order 
behavior 


2 rc^c 


In 


\x - (rp ~ r H ) 
rc 


(5.14) 


Thus the condition (4.8b) is satisfied if b = rc and 


Xq 2) = a ( x + r H) + r c (x + r H ) J 


dx 


(x + r H ) 2 f 


(5.15) 


Substituting these solutions into Eq. (4.12) we obtain 


= 

rH 


rc 

m = 

0 , 

'io = 

ar H 

9 = 

rc 


r H 


(5.16a) 

(5.16b) 

(5.16c) 

(5.16d) 


Since 8$ = 0 and = 1, the Wronskian condition (4.13) is satisfied. Substituting (5.16) 


into (4.15) and (4.17) gives 


Rj = r ° T \ + O(w) , R h = T \ f + O(cj) ,T I = T H = 2 J c ^ h 2 + O(w) (5.17) 

r C + r H r C + r H r C + r H 


and 


r = 


4 r C r H 

( r c + 4) 2 


(5.18) 


in agreement with the results of EH 
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In Sec. IV it was shown that a necessary condition for the gray-body factor to approach 
a nonzero constant in the zero frequency limit is that V e g should change sign somewhere 
between —oo<x*<oo. However, this is not a sufficient condition as we show next by 
considering a massless scalar held with coupling £ ^ 0 to the scalar curvature. It was shown 
in fl2 j that for the l = 0 mode in Schwarzschild-de Sitter spacetime T ~ uj 2 . Using the 
fact that in Schwarzschild-de Sitter spacetime R = 4A, the Klein-Gordon equation takes 
the form 


1 


V=~9 


- £4A</> = 0 


(5.19) 


Transforming to x = 4>/(x + rn) one finds that the effective potential is U e ff = /( ( x +r H ) + 
4£A). For small enough values of |£| its qualitative behavior is the same as in Fig. 4 (see 
Fig. 1 of (I2j). 

Because V e g approaches zero when £ ^ 0 at the horizons in essentially the same way 


as it does when £ = 0, the results in Sec. IV and Appendix A hold. For £ ^ 0, the zero 


frequency mode equation (see (2.6)) is 


d_ 

dx 


d 


(* + r H) f-TZ 


Xo 


dx V x + rn 


- 4£A(x + r H )x o = 0. 


(5.20) 


The crucial feature one loses when considering £ ^ 0 is the existence of a solution for 
Xq 1 ^ which is bounded at x* = Too and which, in terms of the original field variable 

cj) = x/( x + v h), corresponds to a constant field configuration pOJ CEHJ. This can be seen 

_ _ 


for small values of £ by substituting xo = Xo ; = Xo,o + £,\o,i into (5.20). Then to zeroth 


order in £ one finds the general solution (4.18). If Xq 1 ^ is to be regular at both horizons 


then xo,o an d Xo,i must both be regular at both horizons. Examination of (4.18) shows 
that this only occurs for %o,o if b = 0 so that \'o,o = o(x + i"h)- To first order in £ the 
equation is then 


d_ 

dx 


2nd 


(■x + r H ) f— 


Xo,i 


dx V x + rn 


= 4£Aa(x + rn) 2 ■ 


(5.21) 
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The general solution to (5.21) is 


4£Aa 

Xo,i = —— [x + r H ) 


dx' 


(x' + r H ) (*i + r H y 


(5.22) 


f(x') (x' + r H ) 2 f{x') 

with x'] and X2 arbitrary constants. Evaluating the integral near the event horizon (x = 0) 
one finds that the solutions are divergent unless x\ = 0. Evaluating near the cosmological 
horizon (x = rc — Tjj) one finds that the solutions are divergent unless x\ = rc — ru • Thus 
for £ y 0 there is no solution that is regular at both horizons. 


C. BEC acoustic black holes 


This case has similarities and differences with each of the previous cases. As shown 
in m, for the specific model we consider, in which c = c{x) and n, v are constant, the 

II 


metric coefficients in (2.2) are given by / = 


and C± = ^ where as stated in Sec. 


we are ignoring the irrelevant constants m and n. The condensate moves with a constant 
velocity v in the negative x direction and it is assumed that c > |u| for x > 0 and c < |u| for 
x < 0 so there is just one horizon at x = 0 where c(0) = |u|. As for Schwarzschild spacetime, 
the acoustic metric is asymptotically flat since c approaches the constant C 2 > |u| in the 
limit x —»• oo. The crucial difference, however, is in the behavior of the conformal factor 
for the transverse part of the line element C±. In Schwarzschild spacetime Cj_ = x + rn 
is regular at the horizon but diverges at infinity; here Cj_ = -4= is regular both at the 
horizon and at infinity. Thus, as in Schwarzschild-de Sitter spacetime the graybody factor 
will approach a nonzero constant in the zero frequency limit. The effective potential is 


c /-, v s2 " 1,., N ,2 

V efi — ~ 2 1 _ C? C 4 1 _ 1 C 

It is plotted in Fig. 5 for the profile considered in m-- 


(5.23) 


c=\4 + -(4-cj) 


, 2 , (x + d 

1 + - tan" 1 - 

vr V cr„ 


(5.24) 
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FIG. 5. Effective potential for a ID BEC acoustic black hole with the sound speed profile (5.241 
which was considered in [T7] , with (7„ = 8, |t>| = |, c\=\, and = 1. 


with 


d = a v tan 


7r 


Co - c? V u 2 


1 


v 0 ~ o( c l + c 2) 


(5.25) 

L u 2 — C 1 

There is a well with V e Q < 0 located close to the horizon followed by a barrier. The well is 
a generic feature, since near the horizon 


V e ff ~ — (1-2")c < 0 . 


(5.26) 


Thus as expected from the discussion in |I V[ V e s changes sign in the region outside of the 
horizon. 

We next proceed to determine the low-frequency scattering coefficients and gray-body 


factor in this background. The general solution to the zero frequency mode equation (4.18) 
is 


a b [ x '_ c 2 (x) 

Xo = — + — J dx- 


Vc ■ 


c 2 (x ) — v 2 


(5.27) 


The Xq 1 ^ solution is constructed by noting that as x* —> oo, the second term on the 


right in (5.27) has the leading order behavior b^/c^x*. Thus for the condition (4.8a) to be 
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satisfied a = s /02 and 6 = 0 which gives 


(1) / C 2 

Xo = ' — 


c(x) 


(5.28) 


Unlike (5.3), but similar to (5.12), this solution is finite everywhere and approaches a 
constant, -J C 2 /|u|, at the horizon. The condition (4.8b) for x is satisfied if 6 = 1 /y/ty in 
which case 

r „2 


a 


1 


X “ 2> + ^dx)l dX <?-v* 


This solution diverges as \!^ x* at the horizon. Substituting these solutions into Eq. (|4.12h 

y C2 

we obtain 


(5.29) 


^ = \fr\ ’ 

V M 

(5.30a) 

SB = 0 , 

(5.30b) 

if = a , , 

(5.30c) 

9=M. 

V c 2 

(5.30d) 


Since SB = 0 and = 1, the Wronskian condition (4.13) is satisfied. Substituting (5.30) 


into ( |4.15[ ) and ( |4.17| ) gives 

Ri = M 7 C2 + 0(u) , Rh = + o{u) ,T i = T h = + 0(w) . (5.3!) 


M + C2 


M + C2 


C2 + M 


and 


4|u|c 2 
(C 2 + \v\) 2 


(5.32) 


in agreement with the result in [14]. Thus T approaches a constant nonzero value as cu —> 0. 

Examination of the near-horizon dip of the effective potential in Fig. 5 shows that it is 
possible that ‘bound states’ could exist. Such states would come from spatially normalized 
solutions characterized by a purely imaginary frequency uj = iA with A > 0, and they would 
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correspond to classical instabilities (x ~ e A< ) which grow exponentially with time. We have 


found no solutions of this form for the effective potential (5.23) with the profile (5.24). 

To complete this subsection , we point out that for more general models in which c, v 


and n vary with x , we get the metric (2.2) in terms of the spatial coordinate x defined by 
dx = l^-dx, and / = ~ u \ C± = For ID acoustic black holes c, v and n 

(and thus / and C±) always approach finite constants both at infinity and at the horizon. 


From the general solution to the uj = 0 mode equation (4.18) we see that, since Cj_ is 
asymptotically finite, the solution xl;l ^ Fas b = 0. This implies, from (4.19), that B = 0. 
Thus even for this general case, the gray-body factor approaches a nonzero constant in the 
zero frequency limit. 


VI. INFRARED BEHAVIORS OF CORRELATION FUNCTIONS 


The method of finding solutions we have developed in Secs. [TT] - |IV| also allows for a 
detailed analysis of the infrared behavior of correlation functions for a BEC acoustic black 
hole with one spatial dimension HZ! and for a spherically symmetric black hole. In the 
latter case we consider only the s-wave (l = 0) contributions to the correlation functions. 

The gray-body factor enters in the correlation functions through the scattering coeffi¬ 
cients. In the region outside the event horizon, the symmetrized two-point function for the 
quantized massless, minimally coupled scalar held <f> in the Unruh state [16] (the state that 
describes Hawking’s thermal emission) takes the form m 


where 




/ = N 2 / duj± 


xS i x ) Xu* <y)e ^Tc.c. 


J = N 2 du 


sinh(^) 

xi(x)xL*(x')e- Mt - t ' ) +c.c: 


( 6 . 1 ) 


(6.2a) 

(6.2b) 


By making the substitution ^ —> p, where p is the fluid density (see for instance [20]), this applies for 
other systems as well, and not just for BECs. 
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Here j , j stands for the anticommutator and the normalization is N = —k= for a ID 

L J V 4n 

BEC and N = ^ for a spherically symmetric black hole. The mode functions xS and xL 


are defined in Eqs. (3.3). Their asymptotic behaviors are given in Eqs. (3.7), (3.8), (3.10) 


and (3.11). 


We wish to study the infrared behaviors of the integrands of / and J. To do so we 
first use Eqs. (3.7) and (3.10) to write xL and Xu in terms of xY and xT■ The resulting 


expressions for the integrands are 

1 


Tint — 


uj sinh(-7ru;//v) 

i 

C V 

+ 


\T h (uj)\ z x?{x)(x?W)Ye 




(6.3a) 


■4.. = = - {[xTW + fl/MxTW] [(x?V)) - + JJJMtxTM) - ] 

u y. 

l(x?(x)y + RKu)(x?(x ))*] [x?V) + Ri{u)x?{x')\ e^-^} . ( 6 . 3 b) 


The operations of integration over the frequency and taking the limits x* —> ±oo do not 
commute. Thus the integrals must be computed at finite values of x* and only then can 
the limits x* —> ±oo be taken. For fixed x* and x'* in the limit oj —> 0 one has a;lx*I -C 1 


and a;|x M | -C 1. Thus we can use Eqs. (4.11) to write the integrands in terms of %q 
T he result is 


( 1 . 2 ) 


Tint — 


2 K 


7TCU Z 


\ t h(w)\ 2 [x^Wx^V) + 0{u 2 ) 


Tint = - {%o 1) (a;)Xci 1) ( a;, )|l + T2/MI 2 + 0(u ?)} . 


(6.4a) 

(6.4b) 


For Schwarzschild and Reissner-Nordstrom spacetimes, the low-frequency results for the 


scattering coefficients (5.6) show that Tint ~ 0(1) and Ji n t ~ O(w) in the low frequency 


limit. Thus all sources of infrared divergences have been eliminated in (6.2) including the 
factor of - coming from the normalization of the modes and the additional factor of — 


coming from the Planckian distribution in (6.2a). 
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Radically different is the situation for both Schwarzschild-de Sitter and BEC acoustic 
black holes, where (5.17) and ( |5.31[ ) show that there are infrared divergences in I- mt and 
J; n t, i.e. /i nt ~ J 2 and J- mt ~ Q One might be concerned that these infrared divergences 
will persist in other types of correlation functions, such as the point-split stress-energy 
tensor in the gravitational case 
1 


(' T a b{x,t;x ,t )) = - 


4 L 


{g^a'da'd„ + g va 'd a 'd^ - g^g af) 'd a p ({<f>(t, x ), </>(i', x')}) , (6.5) 

where g& are the bivectors of parallel transport (see e. g. [21]) and the density-density 
correlation function in the BEC case (see e. g. M) 

= A(d t d t ' - vd t d x > - vd t >d x + v 2 d x d x >') {{(j)(t,x),(j)(t',x')}) , ( 6 . 6 ) 

where A ~ . Given the structure of the mode functions, each time derivative will 

always bring down one factor of ui, but this is not true for the spatial derivatives due to the 


factors of in (6.1) which give non vanishing u independent contributions. Therefore, 
one would expect that in (6.5) and (6.6) the terms containing two spatial derivatives will 


have the same infrared structure as (6.1) which, in both Schwarzschild-de Sitter and BEC 
acoustic black holes, has a leading infrared divergent term j However, as we show next 
all of the IR divergent terms in (6.5) and (6.6) cancel leaving these correlation functions 
infrared finiteJ3 

For large values of uo one has to be careful because the operations of differentiation 


of (6.1) and integration over u> do not necessarily commute m- However, there is no such 
problem for small values of u. To see what happens to the infrared divergences when spatial 


derivatives of (6.1) are taken we first substitute the solutions for Xg 1 ^ in Scharzschild-de 


Note that in the Schwarzschild de Sitter and BEC cases, by combining the first of (3.91 with (4.6 1 , we 


have Ri = const+ i*0(uj), thus the factor of |l + f?/(u;)| 2 entering in (6.4b I has the form const. + O(oj 2 ). 


8 Recall that for a 4D massless minimally coupled scalar field in Schwarzschild-de Sitter spacetime it is the 
^ = 0 mode that fixes the infrared behavior of the correlation function (the gray-body factor for i ^ 0 
goes to 0 as u> —» 0, see e. g. m ). Similarly, in the (quasi) ID BEC case for modes with non vanishing 
transverse momentum the mode function <f> in (2.1l (and hence also \) has a factor of the form e tk±x± . 
The quantity k± acts as a mass in the ID hydrodynamic theory and naturally regulates any infrared 
divergences [22] . 

9 This result was so unexpected that it was missed in EE! where it was incorrectly stated that there are 
infrared divergent terms in the density-density correlation function. 
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Sitter spacetime (|5. 12) and a BEC analog black hole (5.28), into (6.4) and then divide by 


the factor of C±(x)C±(x') in (6.1). Noting that C±(x) = x + rn in Schwarzschild-de Sitter 
and C±(x) = 1 /\fc{ x) for the acoustic metric, we see that in both cases the factors of C \_ 
are canceled leaving a generic infrared structure of the two-point function of the form 


duu 


duu 


{{(j)(t,x),<fi(t',x')}) = Ci / —j + C 2 / — + IR finite terms , 


ur 


UU 


(6.7) 


with Ci and C 2 constants. It is clear that acting on it with the differential operators in 


(6.5) and (6.6) will remove these IR divergences leaving an infrared finite result (in the 
BEC case this has been checked numerically 


It is also possible using (6.3) to see that the infrared divergences in (6.1) do not con¬ 


tribute to divergences in (6.5) and (6.6) near the horizon by examining the behaviors of the 


modes and xT near the horizon for fixed lu so that one cannot make the assumption 


that a;|a;*| -C 1. Because of the definitions (3.4) it is sufficient to consider xZ aR d X, 


From these definitions and Eqs. (6.3) and (6.1) it is easy to see that for every term in the 


integrands of (6.1) there is always a factor of either xZ ( x ) /Cj_ (x) or xZ( x )/C±(x) and a 


second factor of one of these quantities evaluated at x'. Further in Eqs. (6.5) and (6.6) 


there is at most one derivative with respect to x and one derivative with respect to x' in 
each term. Thus in the limit that we approach the horizon it is sufficient to consider an 
expansion of xZ( x )/C±{x) or xZ( x )/^±( x ) to first order in x. 


For the metrics we consider, near the horizon Cj_has the general form (4.20). In Ap¬ 
pendix B it is shown that near the horizon 


xZ = cos ux* 


+ sin l ox 
xZ = cos UJX* 
+ sin lux* 


Au + aix ^^.o + 0(x 2 ) 


B u + a 1 x 
C u + OL\X 


4 k(k 2 + LU 2 ) 
+ luA^ 


Du + a\x 


4 k(k 2 + LU 2 ) 
tvCu (xDu 

4 k(k 2 + LU 2 ) 
kDu + LuCu 


4 k(k 2 + oj 2 ) 


+ 0(x 2 
+ 0(x 2 ) 
+ 0(x 2 ) 


(6.8a) 


(6.8b) 
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with 


«i = 4 k' 


^1(0) 


(6.9) 


c ± (o) • 

Note that the expression for xZ can be obtained from that for xZ via the substitutions 
A u —> C u and B u —> D u . So we focus on xZ/C±- Combining (4.20) and (6.8a) and 


using (|6.9|) one finds for all values of x 

1 


xZ(x) 

C±(x) 


- { (A u cos xx* + Bcj sin xx *) 


C±( 0) 

+{A U sinwx* — B u cos xx*) 


1 - 


^1(0) 


IX 


C±(0) K 2 + IX 2 


x + 0(x 2 ) 


CKO) 


KX 


jx + 0(x 2 ) 


} . ( 6 . 10 ) 


C_l(0) K 2 + X 2 

Thus terms that after differentiation with respect to x survive in the limit x —> 0 are all 


multiplied by at least one factor of x. So there are no infrared divergences in (6.5) and (6.6) 
in the near horizon limit. 

For completeness we display expressions for xS an d xL i n the near horizon region in the 
Schwarzschild-de Sitter and BEC cases for small values of x but with no assumptions about 


the value of w|x*|. For Schwarzschild-de Sitter spacetime C± = x + rn- Using Eqs. (4.6), 
fl5T6l ), fl5T7l ), flO] ), ( fT8| and ( fiHT] ). we finfl 

= — f e iuJX * [1 + 0{xx)} + R h {x = 0)e~ iuJX * [1 + O(x) + O(wx)]) , 
x + vh th 1 1 


xL 


= T/(w = 0)- 


[1 + 0(x) + 0(xx )] . 


( 6 . 11 ) 


X + fH ' Th 

For the BEC acoustic metric C± = l/y/c(x). Using Eqs. (4.6), (5.30), (5.31), (3.4), (3.8) 


and (3.11), we find 

V~cXu = a/H {e iux * [1 + 0{xx)} + VW\Rh(u = 0)e~ iux * [1 + 0(x) + O(o;s)]} , 
VcxL = VHT^x = 0)e~ luJX * [1 + 0(x) + 0(xx)} . 


( 6 . 12 ) 


10 Note that for Schwarzschild-de Sitter spacetime and the BEC acoustic black hole it is possible to show 
there are no terms of 0{x 2 ) and higher in these expressions which are not also multiplied by some positive 
power of uj. This can be done by first fixing the value of x to be small but nonzero and taking the limit 
uj — ^ 0 in which case one finds that > XT lim w -nj xZ — xK- Next one takes the limit uj — > 0 
in (6.101 and compares with the exact expressions for XoK-U/Cj-l®) using Eq. (5.121 for Schwarzschild 
de Sitter and (|5.28|) for the BEC case. 
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VII. CONCLUSIONS 


In this paper we have developed a rigorous method based on Volterra integral equations 


of the second kind (2.16) to determine the low-frequency behaviors of the scattering coef¬ 


ficients and gray-body factor for a massless minimally coupled scalar field satisfying the 


Klein-Gordon equation (2.1) in the static metric (2.2). These are given in (3.9), (3.12) and 


(3.13), in terms of coefficients (4.6) and (4.12) that can be expressed as horizon boundary 


( 12 ). 

values of the real solutions to the zero frequency mode equation x 0 ’ ’ with the asymptotic 


behaviors (4.8). These results are valid for the modes of the scalar field which are either in 


the zero angular momentum or s-wave sector for a spherically symmetric black hole, or in 
the longitudinal sector of a BEC acoustic black hole with effectively one spatial dimension. 


From the general analysis of Sec. IV we see that there are two qualitatively different 
behaviors for the gray-body factor which depend on the value of the coefficient SS in 


(4.12). If SB ^ 0 the gray-body factor vanishes as cj . As shown in Sec. VA this occurs 


for Schwarzschild and Reissner-Nordstrom black holes. If 38 = 0 the gray-body factor 
approaches a nonzero constant as u —> 0. As shown in Secs. VB and VC this occurs 


for Schwarzschild-de Sitter black holes and a ID BEC acoustic black hole. It was shown 


in Secs. IV and VB that a necessary, but not sufficient condition to have SS ^ 0 is for 


the effective potential, 14ff to change sign. It was also shown in Sec. IV that a necessary 


condition to have SB / 0 is that the zero frequency solution Xq 1 ^ be bounded at the event 


horizon. By definition (4.8a), xd ^ is always bounded either at the cosmological horizon or 
infinity for the cases we consider. 

The results for the scattering coefficients and gray-body factor were used to study 
infrared divergences. The constant gray-body factor at small u for both Schwarzschild-de 
Sitter and ID BEC acoustic black holes implies that Hawking radiation is dominated by an 


infinite number of low energy particles as can be seen from Eq. (1.1). As shown in Sec. VI 


it also implies that the two-point function (6.1) is infrared divergent with the general form 


of the divergences given in (6.7). No such divergences exist in Schwarzschild and Reissner- 
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Nordstrom spacetimes. Despite this fact, we proved that, both away from the horizon and 


close to it, the point-split stress-energy tensor (6.5) in Schwarzschild-de Sitter spacetime 


and the density-density correlation function (6.6) for the ID BEC acoustic black hole are 
always infrared finite. 
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Appendix A: Some results used for the computation of A u 


In this appendix we give several proofs relating to the IR behaviors of the coefficients 


in (2.18). Throughout we focus on A w . The generalization of the proofs to the other 
coefficients is straight-forward. 


We begin by writing the integral in (2.18) for A u in the form 


Aj — 1 + I\ + I2 + -^3 
h = ~ [ dy* sin(ujy*)V eS (y*)xZ(y*) 

^ J — OO 

h = ~ f dy* sin(u}y*)V eS (y*)xZ(y*) 

U J -A 

h = ~ [ dy* sm(uy*)V eS (y*)xl{y*) • (Al) 

w J A 


The first proof involves showing that the first integral vanishes in the limit A —> 00 for all 
uj including u> = 0 if we take the limits u —> 0 and A —> 00 in such a way that ujA <C 1. 
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Note that this condition on the limits allows an expansion of the oscillatory functions in 
the integrand of the middle integral in powers of uy* and also allows one to expand V e g in 
the limit of large negative and large positive y* in the first and third integrals respectively. 


Since the first integral in (Al) covers the region near the horizon it is useful to rewrite 


the general form of the Volterra equations (2.16) for xZ an d Xu so that the integration range 


is —oo < y* < x*. This can be done by substituting (2.17) into (2.16) and using (2.18) 


and (3.2). The result is 


xZ = Aj cos lux* + sin ujx* H— 

UJ 


dy* sin[w(x* - y*)]U eff (y*) xZ(v*) , ( A2a ) 


' —OO 


xZ = C w cosux* + Du sin ujx* + - [ dy* sin[uj(x* - y*)]U e ff (y*) xZ(v*) • ( A2b ) 

^ J — oo 


We bound the first integral in (Al) by bounding the integrand so that 

r -A 


\h\ = 


1 


U! 


dy* sin (ujy*)V e g(y*)xZ(y*) 


<~ [ A dy*\V eS (y*)\ \xZ(y*)\ ■ ( A 3) 

^ ./—no 


Using (A2) and iterating we find 


“ 1 r x * fVu -1 

xZ( x *) = A UJ COS UJX* + B u sin ujx* + V] — / dy{ ... / dy* n sin[w(x* - yj)] 

n _l ^ J — oo J — oo 


xKff(yi)---sin[w(y*_i - y*nWett(y*n) (A, coswy* + B^sinujyZ) . (A4) 


Thus 


\xZ(**)\ <|A,| +\B U \ +f;4 r dy *- f"" 1 d yn\ v es{y\)\---\V e g(yZ)\ (|AJ +|S W |) , 

n= 1 U J -°° 

( °° i rx* rx* \ 

1 + E WAt / - / dy* n \V e g(y* 1 )\...\V eS (y* n )\ , 

n =1 H ' UJ J 

(A5) 


= (|AJ + |-Bj)exp ( - f dyl\V eS (yl] 

\ ^ ./ — no 


We want to use this bound in (A3). Since the integrand in (A3) covers the range —A > 


y* > — oo it is possible to make the further bound 


1 


\xZiv*)\ < (IAJ +1^1) exp ( - / dy\\V e g{y\)\ ) . 


—A 


(A6) 
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Then (A3) becomes 

\h\< (14,1 + | B w |) exp (— [ ^ dytlVM)]) - f A dy*\V eS (y*)\ ■ (A7) 

\ ^ J —oo J ^ J —oo 


This shows that I\ vanishes for all u > 0 in the limit A —>■ oo. For the limit uj 0 
with uA < 1 we can take the limits in such a way that wA 1+e = 1, for some e > 0. Then 
a sufficient condition that I\ should vanish for A —> oo is that 


A 1+e 


r— A 


dy*\V eS (y*) |^0. 


(A8) 


To go further we need to find a more explicit expression for V e g near the horizon. This 
can be done by writing it in terms of a power series in x. Expansions for / and C± in powers 
of x have been given in Secs.[Tl]and VI respectively. Using (2.13) along with (2.8), (2.10), 


and (4.20) one finds 


Ueff = 


= j 

i= 1 


(A9) 


with 


«i = 4 k' 


i £1(0) 

Cl(0) 


(A10) 


Using the dehnition (2.7b) along with the expansion (2.8) for / one can show by iteration 


that near the horizon 


x ( x *) = y~] ftfct 


2knx* 


(All) 


k =1 


Thus it is also possible to write 


VeB = Y'V e2jKX ' ’ 

j =1 


(A12) 


with 


7i = atiPi . 


(A13) 
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If we bound each term in this last expansion and substitute into (A8) then it is easy to 


show that this condition is satisfied for any e > 0 so long as the resulting sum converges P] 
We next bound the third integral in (Al) by bounding its integrand so that 

i noo | noo 

|/ 3 | = - / dy*sm(u;y*)V eS (y*)x c u (y*) < ~ / dy*\V eS (y*)\ \xZ(y*)\ ■ (A14) 

U J A U J a 


Using Eq. (2.16) and iterating we find that 


°° / -j \n no o no o 

xZ(x*) = COS (ux*) + Y -7“ / d Vi- / d V*n sinp(x* - yl)} 

n=l W dx* Jy^_ 1 

x VesiVi) sin[u;(y*_ 1 - y*)] V eS {y n ) cos(wy*) . 


(A15) 


Thus 


lx£(®*)l < exp ( ^ / dyt|U eff (y^)| 


(A16) 


We want to use this bound in (A14). Since the integrand in (A14) covers the range 


A < y* < oo it is possible to make the further bound 

\x c u(x*)\ < exp (- I dyl\V eS (y{)\ 
J A 


(A17) 


Substituting into (A14) one finds 


\h\ < exp f— j dyi|U eff (yJ')A - f dy\V eS (y*)\ . (A18) 

Ja J u J a 

This clearly vanishes for all u > 0 in the limit A — > oo. For the limit u —> 0 with wA <C 1, 
as above we can take the limits in such a way that cjA 1+e = 1, for some e > 0. Then a 
sufficient condition that J 3 should vanish for A —> oo is that 


n oo 

A 1+e / dy*\V e s(y*)\ —> 0 . 
Ja 


(A19) 


In this case however there are restrictions on the value of e which relate to the asymptotic 
behaviors of V e g. For Schwarzschild and Reissner-Nordstrom spacetimes and for a BEC 


Even if the sum is an asymptotic series, each term in this series vanishes in the limit A —> oo. 
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acoustic black hole with sound speed profile (5.24) it is not difficult to show that for large x, 
V e fi ~ X~ 3 . In this case any value of e in the range 0 < e < 1 will work. For Schwarzschild- 
de Sitter spacetime it is not hard to show that V e g ~ e~ 2K ° x * with nc the surface gravity 
of the cosmological horizon. In this case the only restriction on e is that it be positive. 


Next we show how to manipulate I 2 to obtain (4.2). Using the mode equation (2.14) 
one finds that 


1 


h = ~ 

U J -A 


A , * . , (d 2 xZ(y*) . 2 2 , 

dy sin (uy ) + u; x u (y ) 


(dy 


(A20) 


Integrating the first term by parts gives 
~sm(uy*) dxZ(y 


h = 


uj 


dy* 


J -A 


1 

uj . 


dy* 


/ *,dx w (y ) 2 ■ / 

-ucos{u)y )—— -hw sin (uy )x u {y ) 


dy* 


Integrating a second time by parts gives 

"sin(u;y*) dx c u (y* 


h = 

One can then take the limit A 
wA <C 1. 


UJ 


dy* 


— cos(ujy*)Xu(y*) 


J -A 


00 if either uj > 0 or if the limit uj 


(A21) 

(A22) 

0 is taken such that 


Finally we argue that &/ = A + O(ur). First note that from Eq. (2.13) it is clear that 
so long as C± is analytic in x and nonvanishing then Veff is as well. Further, x*(x) is also 
analytic except at x = 0. Thus Xu is analytic both in x* and ui since, as can be seen 


from (A15), the integrand for the integral over y* is analytic in y*, y*_i, and uj. This is 


correct even in the limit uj —> 0. Next note from (A15) that xZ is an even function of uj. 


Then from (2.18) it is clear that A u is also an analytic function of uj and an even function 
of uj. Thus = &/ + 0(uj 2 ). 


Appendix B: Near-horizon computations of xZ an( i XL 

In this appendix we use the Volterra equation to derive the exact near-horizon behaviors 


of the solutions xZ an d X.Z- We begin by substituting the expansion (A12) into (A2a) with 
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the result that 




OO I OO 


Xu( x *) = coswx* + Bu sincaec* + ^ ^ ... ^ 

n=1 [h=i jn =i 


7ji-7j n 

0J n 


dy* 


x sin[ca(x* - y£)]... sin[o;(y*_ 1 - y* n )]e 2nKy ^ ...e 23nKy ^ 
X (Aj cos cay* + B u sin cay*)} . 


C3/n-l 


dy* n 


(Bl) 


Examination of (A2b) shows that the corresponding expansion for xfj is obtained with the 


substitutions A u —> C u and B u D u . 

At this point the integrals can all be computed starting with the integral over y* and 
working in order to the integral over y\. The first integral is 


CVn -1 


dy* n sin[w(y*_ 1 - y*)]e 2>iK2/ ™ (A u cos cay* + B u sincay, 
= wexp(2j„«:y*_ 1 ) (aW cos uy* n -i + sin cay; 


TO—1 1 ) 


(B2) 


with 



(2 j n K) 3 + 4(2 j n K)ca 2 
—jn B u -)- 2caA U) 
(2j n ^) 3 + 4(2j n K)ca 2 


(B3a) 

(B3b) 


The next integral is then 



dy * n -1 sin[ca(y;_ 2 


y*n- i)]e 2(j "- 1+J " )K < 


cos cay*_! 


+ B$ sinejy*_ x ) 


^exp(2(j n _i + jn)^y*n- 2 ) 


coscay n _ 2 


-9 + £(, 2) 


sincay n _ 2 


(B4) 


with A^ and B^ obtained by making the substitutions A —> A^ and B—> B^ 

B^ \ Clearly the same form holds for all of the 


in (B3) followed by A u 


Ai 1 ^ and B u 
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integrals and the result is 




OO I OO OO 


Xu = A w cos lux* + B u sin ux* + ££■■■£ lh-lju exp(2«:(ji + ... + j n 

n=l [jl = l jn = 1 


x cos wx* + B^ sinwx*^ | , 


(B5a) 


Xlj = Coj cos ux* + D u sin ux* + 


OO OO OO 


^ ^ Tji ■■■7j n exp(2K(j 1 + ... + jn) 

n= 1 lil = l jn = 1 


X 


x (cos ux* + sin cox'* 


(B5b) 


If we include terms of 0{x ) but not higher then using (All) and (A13) 


+ B^ smux*^j , 

(B6a) 

+ sinwx*') . 

(B6b) 
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